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Abstract. The LHC data on event-by-event harmonic flow coefficients measured in PbPb collisions at 
center-of-mass energy 2.76 TeV per nucleon pair are analyzed and interpreted within the HYDJET-I--I- 
model. To compare the model results with the experimental data the unfolding procedure is employed. 
The essentially dynamical origin of the flow fluctuations in hydro-inspired freeze-out approach has been 
established. It is shown that the simple modification of the model via introducing the distribution over spa¬ 
tial anisotropy parameters permits HYDJET-I—I- to reproduce both elliptic and triangular flow fluctuations 
and related to it eccentricity fluctuations of the initial state at the LHC energy. 


1 Introduction 

Azimuthal anisotropy of multi-particle production in rela¬ 
tivistic heavy ion collisions is a powerful probe of collective 
properties of sub-nuclear matter created at extremely high 
densities and temperatures (see, e.g., recent reviews [TJH] 
and references therein). It is commonly described by the 
Fourier decomposition of the invariant cross section in a 
form; 

_ d^N 

dp^ 2TTpTdpT'dri 

OO 

1 + 2 '^ Vn{pT,V) cos [n{ip , (1) 

n=l J 

where pt is the transverse momentum, ry is the pseudo¬ 
rapidity, if is the azimuthal angle with respect to the re¬ 
action plane ifh, and are the Fourier coefficients. The 
observation of strong elliptic flow, which is the second har¬ 
monic, V 2 , in heavy ion collisions at RHIC was argued 
as one of the main evidence for strongly-interacting par- 
tonic matter (“quark-gluon fluid”) formation [31I31[S]|B] ■ At 
the LHC, a number of interesting measurements involv¬ 
ing momentum and centrality dependencies of second and 
higher-order harmonic coefficients in PbPb collisions at 
= 2.76 TeV have been done by ALICE [ZUHIIS], AT¬ 
LAS [TUlfTTlfnifninillTS] and CMS [reifTTlfTSlIT^I^ Col¬ 
laborations. In particular, the event-by-event (EbyE) dis¬ 
tributions of second, third and fourth harmonics of the 
anisotropic flow have been obtained |12| . Other impor¬ 
tant observations are the azimuthal anisotropy of jet m 
and charmed meson cinni yields in PbPb collisions, and 
elliptic V 2 and triangular V 3 flow of inclusive 
and identified m hadrons in pPb collisions. 


In our previous study m the second and higher-order 
harmonics of inclusive and identified charged hadrons in 
PbPb collisions at -v/snn = 2.76 TeV were analyzed in 
the framework of HYDJET-I--I- model m- It has been 
shown that the cross-talk of elliptic V 2 and triangular ^3 
flow in the model generates both even and odd harmonics 
of higher order. This mechanism is able to reproduce the 
Pt and centrality dependencies of quadrangular flow V 4 , 
and also the basic trends for pentagonal U 5 and hexago¬ 
nal vq flows. Moreover, it reproduces also specific angu¬ 
lar dihadron correlations including the so-called “ridge- 
effect” m- However, here we restricted ourselves to the 
analysis of the event-averaged harmonics Un(pT)- In re¬ 
cent years, the study of anisotropic flow fluctuations has 
attracted much interest because of their direct connection 
with the geometry of the initial state of a relativistic heavy 
ion collision [5ni5^IB51IMll^[551l571IB51l5^l40ll41lH2] . In the 
present paper, therefore, we analyze the event-by-event 
distributions of the flow coefficients in PbPb collisions at 
the LHC within the HYDJET-I--I- model. 

The paper is organized as follows. The flow fluctua¬ 
tions intrinsic to the HYDJET-I-- 1 - are discussed in Sec. [5] 
Here the probability densities of both longitudinal and 
transverse flow components, as well as the flow modu¬ 
lus, obtained at different collision centralities are shown 
to be nicely fitted to Gaussian. The fluctuations can be 
enhanced by the EbyE Gaussian smearing of the spatial 
anisotropy parameters of the model. Section [3] describes 
the unfolding procedure proposed by the ATLAS Collab¬ 
oration to get rid of the non-flow fluctuations. This pro¬ 
cedure is utilized in Sec.|4]in HYDJET-I—I- calculations to 
compare the model results with the experimental data on 
the same footing. The agreement with the data on fluctu- 
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ations of both elliptic and triangular flow is quite good. 
Conclusions are drawn in Sec. [5j 


2 Inherent flow fluctuations and eccentricity 
fluctuations in HYDJET++ model 


Event generator HYDJET- 1 --I- (the successor of HYDJET 
[13]) is the Monte-Carlo model of relativistic heavy ion col¬ 
lisions, which incorporates two independent components: 
the soft hydro-type state with preset freeze-out conditions, 
and the hard state resulting from the in-medium multi- 
parton fragmentation and taking into account jet quench¬ 
ing effect. The details of this model can be found in the 
HYDJET-I--I- manual [33]. Its input parameters have been 
tuned to reproduce the experimental LHC data on vari¬ 
ous physical observables measured in PbPb collisions [351 
m], namely, centrality and pseudorapidity dependence of 
inclusive charged particle multiplicity, transverse momen¬ 
tum spectra and correlation radii in central PbPb 

collisions, momentum and centrality dependencies of el¬ 
liptic and higher-order harmonic coefficients. 

In order to simulate higher azimuthal anisotropy har¬ 
monics, the following simple modification [381145] has been 
implemented in the model. In original HYDJETH—h ver¬ 
sion m the direction and strength of the elliptic flow are 
governed by two parameters. The spatial anisotropy e{b) 
represents the elliptic modulation of the final freeze-out 
hyper-surface at a given impact parameter &, whereas the 
momentum anisotropy 6 (b) deals with the modulation of 
flow velocity profile. Both 6 (b) and e(b) can be treated 
independently for each centrality, or (basic option of the 
model) can be related to each other through the depen¬ 
dence of the elliptic flow coefficient U 2 (e, <5) obtained in the 
hydrodynamical approach [46] : 


V 2 (e, 6 ) oc 


2 ( 6 -e) 

(I-d2)(l_e2) • 


( 2 ) 


Then, due to the proportionality of V 2 (b) to the initial 
ellipticity £o(b) = bj^RA-, where Ra is the nuclear radius, 
the relation between d( 6 ) and e(b) takes the form [29] : 


^/l + 4B(e + B) - 1 
2B 


B = C(1 — e^)e , e = fceo , 


(3) 

where two parameters C and k are independent on cen¬ 
trality and should be obtained from the fit to the data. 

Compared to the former transverse radius of the fire¬ 
ball, which reproduces the elliptic deformation 


Reii(b, <fi) = Rfo 


-\/l -I- e(b) cos2(^ 


(4) 


Here (p is the spatial azimuthal angle of the fluid ele¬ 
ment relatively to the direction of the impact parameter. 
Rfo is the model parameter which determines the scale 
of the fireball transverse size at freeze-out, and the new 
parameter £ 3 ( 6 ) is responsible for the triangular spatial 
anisotropy. The event plane of the triangular flow, if' 3 , is 
randomly oriented with respect to the plane which is 
fixed to zero in the model calculations. This means that 
the elliptic and triangular flows are generated indepen¬ 
dently, in accordance with the experimental observations. 
Higher flow harmonics are not explicitly generated in the 
model, therefore these harmonics are absent if both V 2 and 
V 3 is absent. 

It should be noted that although the azimuthal aniso¬ 
tropy parameters e( 6 ), 6 (h) and £ 3 ( 6 ) are fixed at given 
impact parameter 6 , they define v-o_(h) only after the av¬ 
eraging over many events due to the inherent model fluc¬ 
tuations. The main source of the flow fluctuations in HY- 
DJET-I—h is fluctuations of particle momenta and mul¬ 
tiplicity. Recall, that the momentum-coordinate correla¬ 
tions in HYDJET-I--I- for soft component are governed 
by collective velocities of the fluid elements, and so the 
fluctuations in particle coordinates are reflected in their 
momenta. The fluctuations became stronger as resonance 
decays and (mini-)jet production are taken into account. 
An event distribution over collision impact parameter for 
each centrality class also increases such fluctuations. 

The detailed study of the event-by-event flow fluctua¬ 
tions is the subject of our present investigation. The possi¬ 
ble further modification of HYDJET-I--I- to match the ex¬ 
perimental data on the flow fluctuations would be smear¬ 
ing of all three parameters, e, 6 and £ 3 , at a given b. 

To get some notion for the inherent model fluctua¬ 
tions, we start with HYDJET-I--I- simulations for the sim¬ 
plest case of central PbPb collisions (b = 0) in which all 
azimuthal anisotropy parameters £(b), 6 (b) and £ 3 ( 6 ) are 
equal to zero. Figure [T] shows the probability densities 
both for each component of the flow vector V„ and for 
its modulus 14 = l^nl, n = 2, 3, 4.0 In this “fluctuation- 
only” scenario the probability densities of are well de¬ 
scribed by two-dimensional (2D) Gaussian functions [JT] 
HE]: 


P(Y„) = ^;^exp[-V'^/(2cr2)] ^ (6) 

whereas the probability densities of 14 have the forms of 
one-dimensional (ID) Gaussians 

p(Vn) = ^exp[-H„V(2cr2)] ^ (7) 

which are obtained from Eq.® by integration over the 
azimuthal angle. 


the altered radius of the freeze-out hyper-surface in az¬ 
imuthal plane takes into account triangular deformation 
as well: 

R(b, (p) = Ren(b, p)[l + £3(b) cos[3(75 - <^£ 3 )]] . (5) 


^ Via V„ we denote the flow vector determined on EbyE 
basis, while the standard definition Vn = (cos [n(p — !^^n)]) is 
reserved for the flow components obtained by the averaging 
over all particles in an event and over all events in the data 
sample. 
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These distributions are characterized by a single pa¬ 
rameter (T„ only, which regulates both the modulus mean 
{Vn) and the width av„ = (V^) — ( 14 ,)^ as 


{Vn) — 

= ^2 - ^ an , (9) 

^ = \/ir7 = 0.523. ,10) 

In HYDJET-I--I- the value of this fitting parameter, 
Gaussian width cr„, appears to be unique for all harmon¬ 
ics with a good enough accuracy: cr 2 — era ~ 0-4 Ri 0.013. 
It depends on a number of model parameters, which were 
already fixed. The main regulator of cr„ is the mean multi¬ 
plicity, which determines the variation of an{b) with cen¬ 
trality. 

Recall, that in all Bjorken-like models with cylindrical 
parameterization the azimuthal anisotropy of the freeze- 
out surface transforms into the azimuthal anisotropy of 
particle momentum distribution proportionally to a term 
{pt sinh Yt/T) cos {(j) — <p) [IH] , arising in scalar product of 
4-vectors of particle momentum and flow velocity of the 
fluid element. Here (j) is the azimuthal angle of the fluid ele¬ 
ment, p is the particle azimuth, T is the freeze-out temper¬ 
ature and Yt is the transverse flow rapidity, respectively. 
The pre-factor before the cosine controls the azimuthal 
angle structure of particle spectrum and its inverse char¬ 
acterizes the fluctuation width squared. We have also veri¬ 
fied numerically that at fixed mean multiplicity in selected 
Pt window the width of anipr) is approximately propor¬ 
tional to the factor y^T/ {px sinh , where is the 

maximal transverse flow rapidity. 

The “true” direction of the flow vector Vn in HY- 
DJET-I--I- for any azimuthal harmonic is pre-defined in 
each event. Therefore, Vnh can be calculated as a longi¬ 
tudinal component of the actual flow vector along this 
known direction, and VqT — as a transverse component 
of Vn perpendicular to the longitudinal one. As shown 
in Fig. [TJ the mean values of the flow components are 
zero, (I4l) = (Kt) = 0 , n = 2,3,4, and, therefore, 
V 2 = V 3 = V 4 as it should be in a trivial case of the flow 
absence in a perfectly central collision. 

Figure [H demonstrates the probability densities both 
for each component of the flow vector Vn and for its mod¬ 
ulus Vn = |V„|, n = 2,3,4 in the case of non-zero flow 
vector (Kil) in PbPb collisions at centralities 20 — 25%. 
Instead of distributions m and © we get here 


p(Kl) 


p(Kt) 

p{Vn) 


1 

\/27rcr2(6) 

1 

\/27rcr2(6) 


exp 

exp 


(Kl - (Kl))^ 

2^^(&)J ’ 


Vn 


K? + (Kl)^ ] / K(Kl) 

\ ° 'v <{h) 


( 11 ) 

( 12 ) 


( 13 ) 


where /q is the modified Bessel fun ction of the fi rst kind 
zero order. Both the width avr, = \J (K?) — (14)^ and the 
modulus mean (P„) are controlled by the true value of 
the flow vector (I4 l) and the width an{b), but cannot 
be cast analytically as functions of (Vhl) and an{b). Note 
also that {Vn) is not equal to (Piil) exactly, and the az¬ 
imuthal anisotropy parameters e( 6 ), 5{b) and € 3 ( 6 ) have 
been tuned earlier at a given impact parameter b in such 
a way that the value of (VnL) extracted from the distri¬ 
bution (ED reproduces just the experimentally observed 
value of Vn{b) entering in Eq.([T|). Similar to the non-flow 
case, presented in Fig. [2 the widths cr„ of the longitudinal 
and transverse distributions shown in Fig. [ 2 ] are approx¬ 
imately the same, ,72 — <73 ~ ,74 r; 0 . 02 , but the distri¬ 
butions become broader. Also, the maxima of piVnh) and 
piVn) distributions are shifted towards zero with rising 
harmonic number n, indicating that V 2 > V 3 > V 4 at this 
centrality. Surely, it is interesting to compare our inherent 
model probability densities, obtained without any addi¬ 
tional special parameters for the azimuthal fluctuations, 
with the experimental data. 

We have also considered including of additional “ec¬ 
centricity” fluctuations in HYDJET-|--(- model. The sim¬ 
plest modification for this purpose is to introduce event- 
by-event Gaussian smearing of the spatial anisotropy pa¬ 
rameters e(b) and £ 3 ( 6 ) with the widths proportional to 
its “unsmeared” values. The coefficients of this propor¬ 
tionality are independent on event centrality and tuned to 
fit the data. Both model versions, with and without the 
smearing, are employed for the extraction of the flow fluc¬ 
tuations and comparison with the available experimental 
data. 


3 Unfolding procedure for flow fluctuations 
analysis 


Unfortunately, the direct comparison of the model distri¬ 
butions given by Eqs. (111111311 with the corresponding ex¬ 
perimental data is impossible. The EbyE distributions of 
anisotropic flow harmonics have been obtained by the AT¬ 
LAS Gollaboration [12] for the distribution of the modulus 
of the flow vector by application of the so-called “unfold¬ 
ing procedure”. The goal of the unfolding procedure was 
to extract the “true” flow vector from the observed one 
by excluding the influence of non-flow effects, such as res¬ 
onance decays and jet fragmentation, as well as the finite 
event multiplicity effect. Therefore, in what follows we will 
show our results before and after the unfolding to be ad¬ 
equate. In order to employ the EbyE unfolding procedure 
for simulated events, the analysis method from m was 
utilized. 

— The EbyE distributions of charged particles in PbPb 
collisions at \/snn = 2.76 TeV with px > 0.5 GeV/c 
and \r]\ < 2.5 are used as input distributions. Fourier 
decomposition of the azimuthal distribution is rewrit¬ 
ten as 
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dN 

dip 


OO 


oc 1 -h 2 V cos [n{ip - <1'°'^^)] = 


n—1 


= 1 + 2 ^ (Ki°x" cos nip + sin np) , (14) 


n—1 


where is the magnitude of the observed per-particle 
flow vector and is the event plane angle. 

The single-particle EbyE distributions are constructed 
and used in our unfolding procedure: 


14°bs ^ yjobs = {cosnp) 

= (sinnp) . 


(15) 


The averaging in last two equations of (|T5|) is per¬ 
formed over all hadrons in a single event. As was shown 
in [ 12 ], the distribution obtained after the unfolding 
procedure did not depend on the method applied to 
obtain Thus, the single-particle method can be 

used for our study. 

The response function is constructed using the “two 
sub-event method” (2SE), namely, the charged parti¬ 
cles are divided into two sub-events with rj < 0 and 
77 > 0. The smearing effects are estimated by the dif¬ 
ference of the flow vectors between the two sub-events, 
for which the flow signal cancels. This distribution is 
fitted to the Gaussian with the width ( 52 SE determined 
mainly by the finite multiplicity effect and the non-flow 
contributions. It was shown in [T2j that the response 
function constructed by such a procedure can be ex¬ 
pressed as 


The difference between the ( 52 SE and (T„ arises mostly 
due to the dynamical flow fluctuations. Therefore, the 
EbyE unfolding analysis excludes the effects related to 
( 52 SE, and leaves the genuine flow fluctuations. 

The application of Bayesian unfolding method for the 
anisotropic flow analysis was checked with heavy ion event 
generators in m- It was shown that the restored density 
distributions were able to reproduce the input 14 , distri¬ 
butions. The non-flow effects were estimated by using the 
HIJING event generator m with and without the imple¬ 
mentation of the flow signal, and were found to be of the 
order of statistical errors. In case of the AMPT model |52j . 
which includes both flow and non-flow fluctuations, the 
difference between the “generated” and “unfolded” flow 
harmonics in semiperipheral Au-|-Au collisions at RHIC 
was found to be small for elliptic flow and significantly in¬ 
creasing in the tails for triangular and quadrangular flows. 
Eluctuations originated from the finite multiplicity effect 
can be evaluated under the assumption of Gaussian mul¬ 
tiplicity distribution [50] : 


Sn = 


(^) 


2(7V)L 


1 + 


/ CTN Y 
\{N)) 


(18) 


where (N) and are the mean value and the width of 
the distribution, respectively. 


4 Comparison of HYDJET++ simulations 
with LHC data 


p(Er^|K)cxEr^exp 


(yobs )2 ^ 


21 


2(52 


<52 


Here S = i 52 Se /2 because we use the full-event 
distribution as an input, and ( 52 SE is the width obtained 
from the difference between the EbyE per-particle flow 
vectors of the two sub-events. 

The constructed response function is used to obtain 
the unfolding matrix: 


M'*' 
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iter 
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„iter 
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^iter+l ^ 


jkCf. 
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(17) 


where Aji is the response function between ej = 14 ’'^® 


(“effect”) and Ci = 14 (“cause”). The true 14 distribu¬ 
tion (cause “c”) is obtained from the measured 14 '’^ 
distribution (effect “e”) using an iterative algorithm. 
Then the Bayesian unfolding procedure is performed 
by means of the RooUnfold package m- 


At first we have checked that HYDJET-I--I- reproduces 
well the experimentally measured correlation between the 
event-averaged elliptic and triangular flow coefficients. The 
results of model simulations are plotted onto the ATLAS 
data m in Fig. [3] One can see that the calculations and 
the data agree well within the 7% accuracy limit. 

Then, we consider the anisotropic flow fluctuations. 
^ 2 ^ 0 ^he simulations and analysis were performed for three 
centrality intervals, namely, 5 — 10%, 20 — 25% and 35 — 
40%, and for two settings of the HYDJET-I—|- model: (i) 
without and (ii) with the additional smearing of spatial 
anisotropy parameters e( 6 ) and £3(6), see Sec. H] The re¬ 
sults for 71 ( 14 ) and p(I4) are presented in Fig.[3]and Fig.jSJ 
respectively, and listed in Table I. Both figures indicate 
that the original version of HYDJET-f-|- (without the 
smearing of e(6) and £ 3 ( 6 )) already includes some dynam¬ 
ical fluctuations due to radial flow. Moreover, the mean 
values of (14) and widths of the distributions av„ in the 
default version of HYDJET-I--I- are quite close to the mea¬ 
sured ones in collisions with centralities up to 25% for the 
triangular flow and up to 45% for the elliptic flow, see Ta¬ 
ble I. Although the agreement can be further improved by 
rescaling of (14) to match the data, such approach would 
be completely misleading. Implementation of the unfold¬ 
ing procedure clearly demonstrates in Figs. I4I5I that the 
initial distributions become more narrow. Thus, the in¬ 
trinsic fluctuations appear to be too weak to match the 
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experimental data m- On the other hand, simple modi¬ 
fication of the model via introducing the normal distribu¬ 
tion over the spatial anisotropy parameters allows us to 
reproduce the measured event-by-event fluctuations of el¬ 
liptic and triangular flow, the distribution widths and the 
event-averaged values of (V 2 ) and (V 3 ). The distributions 
p(Vn) obtained with the set of the smeared out parameters 
are broader in the tails compared with the experiment. 
Unfolding makes it narrower. Here the difference between 
the “initial” and “unfolded” spectra are not so dramatic 
although still noticeable in contrast to that at RHIC en¬ 
ergy, calculated in m within the AMPT model. The two 
additional parameters of the model, appeared in this case, 
are the coefficients of proportionality between the Gaus¬ 
sian widths of the distributions p{e) and ^(£ 3 ) and their 
“unsmeared” values. These two coefficients are fixed to fit 
the data on p{V 2 ) and ^(Va), respectively, for only one ar¬ 
bitrary centrality, whereas for other centralities they are 
the same. 

It is worth noting, that such a simple modification of 
the model also increases EbyE fluctuations for higher or¬ 
der harmonics Vn (n > 3), which arise in HYDJET-I--I- due 
to the presence of elliptic V 2 and triangular V 3 flows, and 
its interference. However, significant sensitivity of high 
harmonic values on their extraction methods makes the 
direct comparison of our simulations with the data even 
more tricky than for V 2 and V 3 . For example, the cen¬ 
trality dependence of quadrangular flow V 4 measured by 
event plane and two-particle cumulant methods is signifi¬ 
cantly weaker than that of V 4 measured by Lee-Yang zero 
method due to large non-flow contribution and increase of 
the flow fluctuations in more central events. Since HYD- 
JET-I—I- was tuned to fit the pr -dependence of V 4 {LYZ}, 
it underestimates V 4 extracted by the event plane or two- 
particle cumulant methods in (semi)central collisions |28] . 
We plan to study the event-by-event fluctuations of higher 
order flow harmonics in the future. 

Few important issues should be clarified still. We uti¬ 
lized normal smearing of the parameters e(b),S{b) and 
£ 3 ( 6 ) in the modified version of the event generator. What 
will happen if the parameters are smeared out with respect 
to a non-Gaussian distribution? To check this possibility 
we opted for a uniform distribution of the key param¬ 
eters within the interval ±(Te to ensure the same mean 
and width values. The distribution p{V 2 ) is displayed in 
Fig- m for centrality bin 20 — 25%, where the signal heav¬ 
ily dominates over the fluctuations. Interestingly enough, 
the generated distribution is very close to the ATLAS 
unfolded curve everywhere, but in the I 0 W-V 2 range. Un¬ 
folded HYDJET-I--I- spectrum, however, is narrow. It re¬ 
sembles the implemented rectangular-shaped V 2 -distribu- 
tion with some rounding of the shoulders because of the 
intrinsic model fluctuations. 

Next question is the Gaussian-like behavior of the ob¬ 
tained spectra after the unfolding. ATLAS Collaboration 
reported some deviations from the Gaussians observed for 
p{V 2 ) distributions in peripheral collisions [15]. The last 
centrality bin in ATLAS analysis is 60 — 65%. At this cen¬ 
trality the default version of HYDJET-I--I-, which works 


reasonably well up to 40 — 45% [55], needs further fine 
tuning in line with other semi-phenomenological models. 
Simply, the linear dependence for U 2 (&) becomes too crude 
here. The new tuned values of the parameters are e = 0.14 
and 5 = 0.25 (cf. with e = 0.16 and 6 = 0.38 in the default 
version). It is worth noting that this is the only modifi¬ 
cation of the model, whereas the ratio (e) /de is kept con¬ 
stant for all centralities in question. Figure 0 shows the 
observed and the unfolded distributions of p{V 2 ) obtained 
in HYDJET- 1 --I- for centrality 60—65% in comparison with 
the ATLAS data. One can see that the model calculations 
agree well with the data. The unfolded distribution pro¬ 
vided by HYDJET-I—I- was also fitted to complex Bessel- 
Gauss (BG) product given by Eq. (|T5)) . Results are plotted 
onto the simulated spectra in Fig. [7] as well. At this cen¬ 
trality the BG fit clearly deviates from the data indicating 
that the model possesses some intrinsic fluctuations which 
cause the distortion of initial Gaussians. This interesting 
problem definitely deserves further investigations. 

Finally, our results can be used for the analysis of fluc¬ 
tuations of the initial anisotropy e„. This approach relies 
on the assumption of a linear response of the flow coeffi¬ 
cient Vn to the corresponding initial eccentricity 


where kn is the response coefficient. For both elliptic and 
triangular flow this assumption works very well, as was 
confirmed by hydrodynamic model calculations [551151] . 
The probability distribution of the flow is connected to 
the initial anisotropy distribution via [55] 

fjf 

p{Vn) = ^ p{eu) ■ ( 20 ) 

dVn 

Inserting Eq. (|T5)) into Eq. (1551) we get 

p{e„) = kn p{Vn) ■ ( 21 ) 


Then, in m the Elliptic Power distribution was proposed 
to parametrize the eccentricity distributions 


P{^n) 


20f6'7T, 

TT 


(1 - £ 2 )“+ 1/2 



( 1 - 4 )“ ^dip 

/1 \ 2 q ;- 1-1 

(1 - COSif) 


( 22 ) 


Here the parameter eg is approximately the mean reaction 
plane eccentricity and a describes the eccentricity fluctu¬ 
ations. In [ 55 ] the authors fitted the ATLAS data on the 
elliptic flow to Eq. ( 1551 ) . It is very tempting, therefore, to 
fit the HYDJET-I--I- generated distributions to Eq. ( 1551 ) 
and compare the extracted parameters, a, £q and fc 2 - The 
fitted curves are plotted onto the model calculations (with 
smearing and unfolding procedure) of ^(£ 2 ) distributions 
for two centralities, 20 — 25% and 35 — 40%, in Fig. [ 5 ] 
Extracted fit parameters are compared in Table 2 with 
those obtained in [ 55 ]. The agreement between the two 
sets is good, indicating that HYDJET-I—I- quantitatively 
reproduces the anisotropic flow fluctuations. 
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5 Conclusions 

The phenomenological analysis of the event-by-event dis¬ 
tributions of anisotropic flow harmonics measured in lead- 
lead collisions at the center-of-mass energy 2.76 TeV per 
nucleon pair has been performed within the two-component 
HYDJET-f-|- model. The unfolding procedure was applied 
to examine the simulated events, thus allowing for the di¬ 
rect comparison of model calculations with the experimen¬ 
tal data. To our best knowledge, this is for the first time 
when the model-generated spectra are filtered by means of 
the unfolding procedure and then compared directly with 
the LHC data obtained by the same method. This proce¬ 
dure removes the non-flow effects, originating, e.g., from 
the decays of resonances and fragmentation of jets, as well 
as the finite event multiplicity effect. 

The essentially dynamical origin of the flow fluctua¬ 
tions in hydro-inspired freeze-out approach has been es¬ 
tablished. The effect is traced to the correlation between 
the momenta and coordinates of final particles and the 
velocities of hadronic fluid elements. The simple modifi¬ 
cation of the model via introducing the distribution over 
spatial anisotropy parameters permits HYDJET-I--I- to re¬ 
produce both elliptic and triangular flow fluctuations in 
heavy ion collisions at the LHC energy. In contrast, an at¬ 
tempt to utilize the uniform non-Gaussian smearing with 
the same (e)/cre ratio failed shortly. The unfolding proce¬ 
dure is sensitive, therefore, to the initial distributions of 
the parameters. It should be implemented in any model 
in case of comparison with the unfolded data. 

For the peripheral topologies the model calculations 
deviate from the Bessel-Gaussian fit to Eq. (nni) thus hint¬ 
ing for some intrinsic fluctuations in HYDJET-I —\- which 
cause the distortion of initial Gaussians. This interesting 
problem deserves to be studied in the future. 
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Gentrality 

5 - 10 % 

20 - 25% 

35 - 40% 


(V 2 } 

0.0522 

0.0976 

0.1186 

ATLAS 

av2 

0.0226 

0.0330 

0.0422 


^V2/(V2} 

0.433 

0.338 

0.356 


{V 2 ) 

0.0489 

0.0853 

0.1213 

HYDJET-b-h 

av2 

0.0170 

0.0231 

0.0353 


^V2/{V2) 

0.347 

0.271 

0.291 


{V 2 ) 

0.0460 

0.0823 

0.1163 

HYDJET-b-h 

(TVs 

0.0081 

0.0095 

0.0147 

unfolded 

^V2/{V2) 

0.176 

0.116 

0.126 


(V 2 ) 

0.0580 

0.0999 

0.1309 

HYDJET-b-h 

(TVs 

0.0240 

0.0371 

0.0483 

(smeared) 

^Vs/(R2) 

0.414 

0.371 

0.369 

HYDJET-b-h 

{V 2 ) 

0.0552 

0.0971 

0.1267 

unfolded 

(TVs 

0.0199 

0.0320 

0.0408 

(smeared) 

^V2/{V2) 

0.361 

0.330 

0.322 


Gentrality 

5 - 10% 

20 - 25% 

35 - 40% 


(^^3) 

0.0281 

0.0344 

0.0373 

ATLAS 

crv3 

0.0147 

0.0178 

0.0192 



0.522 

0.518 

0.513 


(V 3 ) 

0.0341 

0.0428 

0.0555 

HYDJET-b-h 

(TVs 

0.0139 

0.0185 

0.0268 


^vJ{Vs) 

0.408 

0.432 

0.482 


(^ 3 ) 

0.0304 

0.0363 

0.0425 

HYDJET-b-h 

(TV3 

0.0059 

0.0067 

0.0142 

unfolded 


0.195 

0.184 

0.334 


(V 3 ) 

0.0350 

0.0438 

0.0527 

HYDJET-b-h 

(TV3 

0.0178 

0.0225 

0.0276 

(smeared) 


0.509 

0.514 

0.524 

HYDJET-b-h 

(V 3 ) 

0.0297 

0.0356 

0.0387 

unfolded 

(TV3 

0.0148 

0.0180 

0.0203 

(smeared) 


0.499 

0.504 

0.526 


Table 1. Mean values (I4i), widths and ratios crv„/(V„) 
(n=2,3) for HYDJET-I--I- simulations and ATLAS data |12 |. 



Gentrality 

20 - 25% 

35 - 40% 


a 

56 ±6 

24 ±3 

ATLAS 

eo 

0.25 ± 0.02 

0.31 ±0.08 

(reanalysis in |55p 

k2 

0.40 ± 0.02 

0.35 ±0.01 


a 

48 ±7 

35 ±3 

HYDJET-b-h 

£0 

0.25 ± 0.02 

0.30 ±0.01 

(smeared and unfolded) 

k2 

0.40 ± 0.02 

0.36 ±0.02 


Table 2. Parameters of the fit of eccentricity distributions to 
Eq. (1221) . See text for details. 
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Fig. 1. The probability density distributions for two components of the flow vector, VnT and Fql, and for its modulus Vn = 
n = 2, 3, 4 in the case of zero signal in HYDJET-|--|- (central PbPb collisions at impact parameter b = 0). Dashed curves 
on the left and middle plots show two-dimensional fit of simulated HYDJETH—h points to Eq. dashed curves on the right 
plots present the ID-fit to Eq. 0. 





Fig. 2. The same as Fig. [T] but for non-zero signal in HYDJET-I--I- (20-25 % centrality of PbPb collisions). Dashed curves in 
the left, middle and right columns indicate the fit of simulated HYDJET-I—I- points to Eqs. (HD), HD and m, respectively. 
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Fig. 3. The correlation between the event-averaged elliptic flow V 2 and triangular flow V 3 of charged hadrons at transverse 
momentum 0.5 < pt < 2 GeV/c and pseudo-rapidity |t;| < 2.5 for ten 5%-centrality intervals in the centrality range 0 — 50% 
of PbPb collisions at v^nn ~ 2.76 TeV. Closed circles denote ATLAS data from [T^, asterisks represent HYDJET-I—I- events. 
Line is drawn to guide the eye. 
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V2 V2 V2 


Fig. 4. The probability density distributions of elliptic flow V 2 in three centrality intervals 5 — 10% (left), 20 — 25% (middle) 
and 35 — 40% (right). Dashed and solid histograms present the results for simulated HYDJET-I—I- events before and after the 
unfolding procedure, respectively. The top/bottom row shows the model results with/without the additional smearing of spatial 
anisotropy parameters. The closed points are ATLAS data from m- 



Fig. 5. The same as Fig. |T]but for the triangular flow V 3 in three centrality intervals. 
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Fig. 6. The probability density distributions of elliptic flow V 2 in centrality interval 20 — 25% obtained with uniformly smeared 
parameters e and S. Dashed and solid histograms present the results for simulated HYDJET-I—I- events before and after the 
unfolding procedure, respectively. ATLAS data are shown by full circles. 



Fig. 7 . The same as upper row of Fig. [4] but for the centrality interval 60 — 65% only. Dashed and solid histograms present the 
results for simulated HYDJET-I—I- events before and after the unfolding procedure, respectively. ATLAS data are shown by full 
circles. Solid curve shows the Bessel-Gauss fit to Eq. m- 
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Fig. 8. The probability densities P{en) obtained in HYDJET-I—I- calculations (with smearing and unfolding) of PbPb collisions 
at y^SNN = 2.76 TeV with centrality 20 — 25% and 35 — 40%. Full lines represent the fits to the Elliptic Power distribution (I22II . 






